Observing different phases for the dynamics of entanglement in an ion trap 
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The evolution of the entanglement between two oscillators coupled to a common thermal en- 
vironment is non-trivial. The long time limit has three qualitatively different behaviors (phases) 
depending on parameters such as the temperature of the bath {Phys. Rev. Lett, 100, 220401). 
The phases include cases with non-vanishing long-term entanglement, others with a final disentan- 
gled state, and situations displaying an infinite sequence of events of disappearance and revival of 
entanglement. We describe an experiment to realize these different scenarios in an ion trap. The 
motional degrees of freedom of two ions are used to simulate the system while the coupling to an 
extra (central) ion, which is continuously laser cooled, is the gateway to a decohering reservoir. The 
scheme proposed allows for the observation and control of motional entanglement dynamics, and is 
an example of a class of simulations of quantum open systems in the non-Markovian regime. 



I. INTRODUCTION 

Entanglement is not only an essential feature of the 
quantum world but also a physical resource enabling the 
manipulation of information in non-classical ways p[]. 
Understanding the evolution of entanglement for open 
systems is thus important from both a fundamental and 
a practical viewpoint. In fact, developing novel quantum 
information technologies requires a thorough characteri- 
zation of the process of decoherence by which entangle- 
ment is typically degraded. 

Considerable efforts were made in this direction and 
it was shown that even simple quantum open systems 
exhibit non-trivial features in the evolution of entangle- 
ment. For instance, it was observed that entanglement 
can vanish in finite time even if coherences only decay 
asymptotically j2j . The long-time behaviour of entangle- 
ment for qubits interacting with bosonic baths has been 
analyzed in [3] , where different asymptotic regimes where 
identified. Recently, a system formed by two resonant 
oscillators coupled to a common thermal reservoir was 
found to exhibit three different dynamical phases char- 
acterizing the asymptotic entanglement 

In this paper we propose and analyze an experiment to 
observe these three phases, using an array of three cold 
ions in a linear trap. This experiment can be viewed as 
a simple instance of a quantum simulation of the evolu- 
tion of quantum open systems in the non-Markovian and 
non-perturbative regime. Trapped ions simulators were 
introduced in [5] and are currently under active investi- 
gation [BVfTU] . Our proposal builds on these ideas in order 
to study the dynamics of entanglement between motional 
degrees of freedom. We note that motional entanglement 
has recently been demonstrated in ion traps in a static 
regime [TT]. 

Let us first describe the problem studied in (4j and then 
show how its essential features could be reproduced in an 
ion trap. The subsystems are two equal harmonic oscil- 
lators identically coupled to a common oscillator bath 
(the coupling is bilinear in position). In the long time 



limit there are three qualitatively different possible be- 
haviors (phases), depending on parameters such as the 
bath temperature, the squeezing of the initial state, etc.: 

1) The system may end up containing a non- vanishing 
amount of entanglement for all sufficiently long times, or 

2) there may be an infinite sequence of events of disap- 
pearance followed by revivals of the entanglement, or 3) 
the asymptotic entanglement may be exactly zero. 

These three possibilities exist in this simple system but 
have never been observed in the laboratory. Recently, a 
scheme has been proposed to observe two of these phases 
in a cavity-QED implementation [12 . Our goal here is to 
propose an experiment in which by varying accesible pa- 
rameters the three different behaviors may be attained, 
and detected. The general idea is illustrated in Fig [ij We 
consider a linear harmonic trap with three ions. The two 
ions at the ends of the chain are assumed to be equal and 
will constitute the subsystems, while the central ion will 
provide an environment with which the subsystems inter- 
act. In turn, the motion of this ion will be laser-cooled 
(pumping energy out of the normal modes in which this 
ion takes part). We will show that these ingredients are 
sufficient for the appearance of the three entanglement 
phases. Either the axial or transverse degrees of freedom 
can be considered; we analyze first the case of transverse 
modes, and mention later the differences that arise if one 
uses axial modes. 

The paper is organized as follows: In section [TTlwe 
introduce our system and its dynamics. Section |III| is 
devoted to the asymptotic state and its entanglement 
properties. In section |IV| we consider a family of initial 
states and find the corresponding phase diagrams for the 
asymptotic entanglement. Section [V] explains how the 
different behaviours can be detected. Finally, we con- 
clude in section IVTl 

II. DESCRIPTION OF THE SYSTEM 

We consider a linear chain of three ions, neglecting 
the motion along one of the transverse directions {y) in 
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FIG. 1. A linear chain of three trapped ions. The ions at the 
ends will be equal and correspond to the subsystems. The 
central ion, coupled to a cooling laser, will be the source of 
decoherence of the system. 



which the trapping is assumed to be tight. The motional 
Hamiltonian, in terms of their longitudinal and trans- 
verse coordinates (z and x respectively), is [131 [T4] : 
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Here, m is the mass of the ions at the ends of the chain, 
while LOz and are the trapping frequencies for these 
ions in the axial and transverse directions respectively. 
The transverse trapping frequency scales inversely with 
the mass of the ion, so it will be different for the central 
ion, which will be of a different species [TS]. We choose 
frequencies so that the stable configuration is a linear 
array, and the harmonic approximation is valid. The 
equilibrium distance between ions is then given by: 
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and the linearized Hamiltonian for the transverse motion 
is: 
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ratio u!x/(^z, the strength of the system-environment cou- 
pling is modified (when Ux 3> Wz the normal modes are 
approximately equal to the local modes so the coupling 
tends to vanish in this limit). 

The use of collective coordinates x± — {xa ± xb)/^/2 
is convenient since the center of mass of the system, x+, 
is coupled to the environment while the relative motion, 
X-, is not. In this way the model is analogous to the one 
studied in [l] , the main difference being that the environ- 
ment consists of a single oscillator. However, the effective 
size of the environment can be enlarged by coupling the 
central ion with a laser tuned to cool down the modes 
in which this ion participates. This produces an effect 
similar to that of a low-temperature reservoir. 

A useful approximation for the effect of the laser cou- 
pling can be obtained in interaction picture and assuming 
that the system is in the Lamb-Dicke regime. Then, if 
the laser is tuned to the red sideband of one of the mo- 
tional modes (n, where n labels either the center of mass 
or the third mode), after a rotating wave approximation 
we find the Hamiltonian |16j : 



H 



(int) 
L-I 



Ml 



2mv, 



-{a 



h.c.) (6) 



Here fifl is the Rabi frequency, k is the component of 
the wave vector along the mode direction, and tp is the 
phase of the laser. a± are the transition operators for 
the internal state of the ion, cr+ = |e)((7|, cr_ = |(?)(e|, 
while a„,ajj are annihilation and creation operators for 
the normal mode. 

Through this coupling, excitations can be transferred 
between the motional modes and the internal state of 
the ion. We shall assume that \g) is a stable internal 
state, while |e) decays back to \g). Thus, the effect of the 
coupling to the laser will be to cool down the motional 
mode n |T7j. We shall assume that each of the modes in 
which the central ion participates is damped in this way. 
The setup we are considering is thus a particular case 
among the ones studied in [18] for sympathetic cooling. 

III. ASYMPTOTIC ENTANGLEMENT 
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Because of the symmetry, the normal modes of this 
Hamiltonian have well defined parity: there are two even 
modes, which involve the motion of all ions, and one odd 
mode in which the central ion is at rest. 

As mentioned before, we associate the two ions at the 
ends of the chain with two susbystems A and B. The cen- 
tral ion provides an effective environment, and the two 
subsystems are equally coupled to it. By modifying the 



As in we shall restrict our analysis to Gaussian 
states. A Gaussian state p is fully characterized by the 
values of the first moments, Tr(pi?j), and the covariance 
matrix: 
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where the curly brackets denote the anticommutator. 
The entanglement properties are only determined by the 
covariance matrix C, since changes in the first moments 
can always be achieved by local operations. 

According to our simple model for laser cooling, the 
odd mode of the chain, corresponding to coordinates a;_, 
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P-, evolves following its own free dynamics. In contrast, 
the even modes approach the vacuum. Thus, the second 
moments of a; 4-, p+ approach equilibrium, characterized 
by {{x+,p+}) = and: 
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Here, el and e2 label the two even modes, ujei/e2 &re 
the corresponding frequencies, and Cei/e2 the amplitude 
coefficients transforming x+ to the even normal modes. 

From the covariance matrix we can calculate the en- 
tanglement between A and B, quantified by the logarith- 
mic negativity En [T^]. If the odd mode is squeezed, 
the entanglement between subsystems may oscillate as 
this mode evolves. The logarithmic negativity reaches 
its maximum and minimum values when the dispersions 
Ax-, Ap_ are extremal (i.e. when the odd mode squeez- 
ing is along the position or momentum axis) and at these 
times it is given by 
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More generally, the asymptotic behaviour of entangle- 
ment is determined by the following three quantities char- 
acterizing the final state [J : 
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Scru = 2 In [j^Ax+Ap+Ax^Ap^j , (14) 

where lUq^^ is the frequency of the odd mode. Here, r 
is the squeezing of the odd mode, Vcrit is related to the 
asymptotic squeezing of the system's center of mass, and 
Scrit to the final entropy of the system formed by the 
ions A and B. These parameters depend on the choice of 
ions and the ratio between the axial and transverse fre- 
quencies, as well as on the initial state of the odd mode. 

The possibility to observe the three dynamical regimes 
in our model is conditioned by the values that these three 
parameters may take [1| . Typically, entanglement will be 
large for large squeezings, and will decrease with increas- 
ing entropy. For instance, the entanglement between sub- 
systems when the chain is in its vacuum state is given by 
the difference: 
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the value for the entropy of the system formed by ions A 
and B when the chain is in its ground state (in general 



Scrit > Smin, for a given choice of masses and frequen- 
cies). In order to obtain large values for rcrit, we need 
the coupling between ions to be strong, which in turn im- 
plies being not too far away from the transition between 
linear and zig-zag equilibrium configurations |il4i . On 
the other hand, the parameters must be chosen in such 
a way that the harmonic approximation is valid (i.e. not 
too close to the transition). To achieve good values for 
E'jtj , it turns out convenient to have a central ion lighter 
than the other two. By using different ion species for the 
subsystems and the environment one also reduces the un- 
wanted coupling of the laser to the subsystem ions. For 
the case in which A, B are ^^Mg ions and the central ion 
is ^Be, Fig. [2] displays the values of rcrit, Smin and i?^ 
as functions of uir^/uj^- We note that these species have 
been used to demonstrate motional entanglement in |llj . 




FIG. 2. Values of the parameters E% (red, continuous), rcrit 
(black, dashed), and Smin (blue, dash-dotted) as functions 
of the ratio lox/l^z for a chain of two ^*Mg ions with one 
^Be ion placed between them. For values of uj^juiz smaller 
than 1.673. . . , as indicated by the vertical line, the linear 
configuration becomes unstable. 



IV. STATE PREPARATION AND 
ENTANGLEMENT PHASES 

In order to make the different asymptotic situations 
actually observable one should find a procedure to pre- 
pare a family of suitable initial states. Conceptually the 
simplest initial state is the vacuum of the chain Hamil- 
tonian which, as pointed out in [20], is a squeezed en- 
tangled state in terms of the local oscillators A and B. 
In this case, the system does not evolve in time, and the 
entanglement is given by i?^. More generally, one can 
start with a thermal state instead of the vacuum. This 
also corresponds to a Gaussian state and leads to a sit- 
uation with a larger final entropy, Scrit > Smim without 
changing the values of r and rcrit- Actually, the asymp- 
totic state depends on the initial temperature for the odd 
mode only. By varying this parameter, associated to the 
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mean population of the mode, the final entropy is 
modified. Thermal states were experimentally obtained 
in |21j from Doppler cooling of a single trapped ion. 

The squeezing of the initial state can in turn be var- 
ied by an abrupt relaxation of the transverse trapping 
frequency, uj'^ ^ oj^ = ^'xl f- After this operation the 
dispersions remain untouched, and since the frequency 
^odd changes, the squeezing of the odd mode is modi- 
fied to r = ^ln(/). Scrit and Vcrit depend on the final 
value of ujx but not on the expansion factor /. More 
sophisticated ideas concerning the creation of squeezing 
using time modulation of the trapping potentials were 
discussed in [22j. An alternative method was presented 
in |21j . where the squeezing in the motion of a single 
trapped ion was experimentally induced obtaining values 
r w 2. This was achieved by the coupling to a laser with 
a detuning equal to twice the frequency of the mode, and 
the same procedure could be used here to squeeze the 
odd mode. 

One can in this way modify the two relevant parame- 
ters Scrit and r, while the value oircnt is only determined 
by the transformation to the normal modes. It is thus 
possible to sample over the three asymptotic regimes, as 
shown in Fig. [s] (for the case lOx = 1.7 lOz)- The be- 
haviour of the asymptotic entanglement as a function of 
time is illustrated in Fig. |4]for different initial states. 
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FIG. 4. Different asymptotic behaviours for entanglement 
(measured by the logarithmic negativity) as a function of 
time. The choice of ions and frequencies is the same as in 
the previous figure. The cases plotted correspond to differ- 
ent values of the parameters {n_) for the initial population 
of the odd mode and r, the subsequently induced squeezing: 
(n_) = 0, r = (black, continuous); (n_) = 0, r = 0.3 (red, 
dashed); (n_) = 0, r = 1 (green, dotted) and {n_) = 1, 
r — 0.7 (blue, dash-dotted). Time is in units of lOz^ , and the 
time origin is arbitrary. We note that the entanglement of a 
Bell pair would correspond to En — 1- 
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FIG. 3. Asymptotic behaviour for the entanglement in a 
Gaussian state parameterized by the initial mean population 
of the odd mode, {n_), and the subsequently induced squeez- 
ing of the mode, r. The width of the area corresponding to 
the phase with death and revivals is determined by the value 
of rcrit, which depends on the choice of ions and on the ratio 
uJx/uiz- The case plotted corresponds to a chain of three ions, 
one ^Be ion in the middle and two ^^Mg ions at the ends; the 
trapping frequencies are chosen in the form lOr^ = 1.7 uj^- 



The axial motion of trapped ions could be used with 
similar results, though the typical values for Tcrit (and 
accordingly E^) are smaller. For our choice of ions, the 
relevant parameters for the axial case are given by Tcrit — 
0.14, Srain - 0.04, Compared to Vcnt - 0.30, Smtn - 
0.12 for the transverse motion with ojx = 1.7 lo^ as in 
Figs. [3] and [4j When the axial modes are considered. 



changing the value of the mass ratio between ions does 
not alter the results much: the phase diagram for our 
choice of ions, shown in Fig. [5] is very similar to the one 
obtained for three equal ions. We note that in this case 
the technique of |5T] seems to be the right strategy to 
generate squeezing, as the sudden alteration of trapping 
frequency would modify the equilibrium positions in such 
a way that the harmonic approximation could only hold 
for very small values of squeezing. 



V. DETECTION 

For the experimental observation of the three entan- 
glement phases, it is enough to measure the asymptotic 
dispersions Aa;+, Ap_|_, and the time- varying quantities 
Ax_, Ap_. The discrimination of the different phases 
only requires the evaluation of the negativity (11) for the 
times when the dispersions take extremal values. The 
necessary variances can be inferred from AQn: AP„, with 
Qn = {an + ajj)/v^, Pn = i{an - a\^_ 
less coordinates for the normal modes. 



ai^) /\/2 the dimension- 



From blue-sideband oscillations together with fluores- 
cence measurements, it is possible to reconstruct the pop- 
ulation of the different levels for each mode [15]. Thus, 
when the first moments vanish the covariance matrix C 
for each mode can be determined in the following way: 
Firstly, the mean phonon number is related to the trace 
in the form: 



Tr(C) = 2{a^a) + 1 



(17) 
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FIG. 5. Asymptotic behaviour for the entanglement in a 
Gaussian state parameterized by the initial mean population 
of the odd mode, and the subsequently induced squeez- 

ing of the mode, r, when the axial modes are considered. The 
case plotted corresponds to a chain of three ions, one ^Be ion 
in the middle and two ^^Mg ions at the ends. 

Then, let us assume that we apply a laser pulse inducing 
a force which changes momentum in a value 6P. The 
probability to find the mode in the ground state after- 
wards is given by 

F{SP} = Vdet(M) e-(^^)'^'^^^/^ (18) 

with 

= ^ + C- (19) 



From this set of measurements one can find C (up to an 
unimportant sign of the off-diagonal elements) for each 
mode at each given time. In this protocol, only the force 
needs to be faster than the chain dynamics, since the 
populations are not modified by the free evolution. We 
note that the ability to perform measurements of time- 
dependent motional observables has recently been de- 
mostrated in [24 , allowing for the quantum simulation 
of the Dirac equation. 

VI. CONCLUDING REMARKS 

In conclusion, we have shown how the symmetry and 
interactions in the three-ion chain allow us to reproduce 
the essential features of the system studied in includ- 
ing the three different entanglement phases. Further- 
more, the observation of the asymptotic entanglement 
dynamics seems possible with available ion-trap tech- 
niques. As a last remark, we note that a more detailed 
investigation of the dynamics of the laser cooling process 
would allow for the description of the approach to the 
asymptotic regime, in terms of a non-Mar kovian master 
equation. The non-Markovian character of the process 
can be varied depending on the effective temperature of 
the bath and strength of the coupling between ions. The 
latter can be modified by the choice of trapping poten- 
tial, while the former depends on the laser detuning. This 
would allow for the extension of the ideas in [25 for the 
case of more than one particle, and will be the subject of 
future research. 
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